Critical disorder effects in Josephson-coupled quasi-one-dimensional superconductors 
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Effects of non-magnetic randomness on the critical temperature T c and diamagnetism are studied 
in a class of quasi-one dimensional superconductors. The energy of Josephson-coupling between 
wires is considered to be random, which is typical for dirty organic superconductors. We show that 
this randomness destroys phase coherence between the wires and T c vanishes discontinuously when 
the randomness reaches a critical value. The parallel and transverse components of the penetration 
depth are found to diverge at different critical temperatures an d T c , which correspond to 

pair-breaking and phase-coherence breaking. The interplay between disorder and quantum phase 
fluctuations results in quantum critical behavior at T — 0, manifesting itself as a superconducting- 
normal metal phase transition of first-order at a critical disorder strength. 

PACS numbers: 74.78.-w, 74.62.-c, 74.70.Kn, 74.50. +r 
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Quasi-onc-dimensional (quasi-lD) organic conduc- 
tors, including the charge-transfer (Bechgaard) salts of 
(TMTSF) 2 X (where TMTSF stands for tetramethyl- 
tetraselenofulvalinium and X = PFq,CIOa, NO3 being 
a strong electron acceptor or anion) [1] and A-15 com- 
pounds [2], attract enhanced interest since the discovery 
of superconductivity in {TMTSF^PFq. Low tempera- 
ture properties of the organic superconductors are very 
sensitive to disorder. Alloying anions, x-vay irradiation 
and cooling rate controlled anion reorientation introduce 
non-magnetic randomness into the system while leaving 
the backbone structure and the unit cell of the organic 
superconductors to a large extent unchanged. There is a 
common agreement that disorder, introduced by means 
of these experimental methods, must be characterized as 
non-magnetic, and yet it was shown [3-5] to suppress the 
superconducting (SC) phase. 

Effect of disorder on the SC phase has a long-standing 
history. According to the Anderson's theorem [6], the 
SC critical temperature T c for s-wave pairing is insen- 
sitive to the scattering rate on non-magnetic impurities. 
Magnetic impurities break time-reversal symmetry of the 
s-pairing, suppress at the same time the SC phase [7]. 
Strong disorder of non-magnetic impurities may how- 
ever destroy (i-wave pairing [8]. Interplay between su- 
perconductivity and Anderson localization in a strongly 
disordered superconductor was shown [8-17] to result in 
spatial inhomogeneity of the order parameter. High pu- 
rity of the organic superconductor backbone even in the 
dirty limit seems to exclude a spatial inhomogeneity of 
the order parameter modulus along SC wires, offering 
an opportunity for another mechanism of disorder-driven 
superconductor-normal metal phase transition. Effects of 
order parameter phase fluctuations on T c have also been 
studied in low-dimensional superconductors [18-21]. It is 
well known that there is no SC phase transition in ID and 
two-dimensional (2D) systems [22], since strong fluctua- 



tions of the order parameter phase destroy off-diagonal 
long-range order (ODLRO) in a single SC wire or film. 
Strong phase fluctuations in clean quasi- ID supercon- 
ductors have been shown [18, 19] to suppress T c below 
a mean-field transition temperature. Classifying the su- 
perconductors with small stiffness as bad metals, Emery 
and Kivelson have evaluated [20] a critical temperature 
of phase ordering by formally dividing a clean bulk su- 
perconductor into small regions with well defined phase, 
and have shown strong suppression of SC phase by phase 
fluctuations. Nevertheless, effects of disorder on phase 
fluctuations are neglected in all of these papers. 

In contrast to these previous activities we study in 
this Letter a suppression of superconductivity as a result 
of the destruction of the order parameter phase coher- 
ence by disorder. We consider weakly linked quasi-one- 
dimensional superconductors with random Joscphson- 
couplings between pure one-dimensional (ID) SC wires. 
Singlet pairing is considered to occur inside the wires. 
Therefore, we assume that non-magnetic randomness 
does not affect the order parameter amplitude. We 
demonstrate in this Letter that (i) non-magnetic ran- 
domness in the Josephson-coupling destroys correlation 
of the phases between different chains in quasi-lD super- 
conductors even in the classical phase fluctuation regime, 
(ii) randomness yields quantum critical behavior. A SC- 
normal metal phase transition occurs at T — with in- 
creasing the strength of disorder, and that (iii) a suppres- 
sion of the SC phase occurs discontinuously in both classi- 
cal and quantum phase fluctuations as a first-order phase 
transition when the disorder-strength reaches a critical 
value. We derive that parallel and perpendicular com- 
ponents of the penetration depth, Ay and A_l diverge at 

different critical temperatures T c and T c , which corre- 
spond to pair-breaking in the wires and to phase coher- 
ence breaking between the SC wires, respectively. 

Classical fluctuations of the phase. The free energy 
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functional of a quasi-lD superconductor weakly linked 
with Josephson coupling energy -Ejj+g between nearest- 
neighbor chains can be written in the presence of the 
magnetic field B as 
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where <p$(z) denotes the order parameter phase, A = 
{A^,A Z } is the vector-potential, and ivi^T) — 
iVi 1 -*(0)r(T) is the linear density of SC electrons with 



t(T) = and 7V S (1) (0) = ~ ^ ^ ^ 



If at T < T r l 



Dimensionless coordinates r = {j,z} are introduced on 

the scale of longitudinal £11 = h Ns and transverse 

£_l ~ components of the coherence length. We assume 
the Josephson energy -Ey+g to be a random parameter 
with Gaussian distribution given by 



P{E i>i+s } 



^expj- (£?JJ+g ~ £?g)2 j. (2) 
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Employing the replica trick one can integrate out the 
Gaussian disorder to obtain the average value of the free 
energy T = — T(lnZ) as 

t = -t / ft — =^ ( /s; i s - 
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with F = iVj 1 ^ E / dz 
j 

£> g - ^ (1) Cn^Cj, g )[i - cos(^ - <p J+g )]\. 
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where £j jg is a Hubbard-Stratonovich auxiliary field. The 
average value of a given functional C({<pj}), e.g. cos<pj 
or cos(<£>j — <Pj+ g ), can be obtained according to the re- 
lation «C({¥>j})» - -T^ (In Z)\ m=0 by adding the 

source term JZj / ^'^({Vj}) to the free energy func- 
tional, which yields for the correlator 
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where the double bracket ((...)) means averaging over 
thermodynamic fluctuations and over randomness. In 



order to estimate an asymptotic behavior of the cor- 
relator, e.g. ((cos(</9j — ifj+g))) we write the in- 
tegrand of Eq.(5) as exp{— A^i 1 ^ 2 ^ 2 /(Cj lg )}, and ap- 
ply the stationary-phase approximation to determine 
an extremal value of the auxiliary field (j ig minimiz- 
ing the function /(Cj.g)- The saddle point value of 
Cj,g is obtained to be Q. s = / dz((cos(ipj(z) — 



constant ivi 1 ^! on the exponent can be estimated to be 
-^j ~ 10 3 for the organic supercon- 
ductors with ep being the Fermi energy, which ensures a 
sharply peaked saddle point of the integrand. The ther- 
modynamic averages in the expression of (j ig are taken 
with the free energy functional, given by Eq.(4) at the 
saddle point (j.g = (j,g- So, a contribution of the non- 
magnetic randomness to the effective free-energy func- 
tional is proportional to the variance of the phase cor- 
relator, which gives an idea on the form of the disorder- 
dependent term in the effective functional. 

The critical temperature for the quasi-lD supercon- 
ductors can now be found from Eq. (5), written for cost^j 
by using the self-consistent mean-field method [18], which 
consists in replacing the phase correlations of the co- 
sine term by J2 S E %\ 1 ~ cos(^O) - f i+s (z))} -> E ± [l - 
({cos(tp))) eff cos(<p(z))], where E ± = J2 s E g- For a 
clean system ((cos(<p))) e // was chosen [18] to be equal to 
(cos(y)). For the disordered superconductor we choose 

«cos^» e// = ((cosy))-^, ^-^-^ 08 ^ 2 . Such 
a form of ((cos <p)) e ff is similar to the expression of 
the saddle point value for the averaged order param- 
eter. Taking advantage of the smallness of (E± — 
iVj 1 ^! WQ((cos(ip))) e ff near T c , we expand both the nu- 
merator and the denominator of the integrand of Eq.(5), 
written for ((cos(y))) e //, in this parameter. The thermo- 
dynamic averages become pure one-dimensional after this 
expansion, which can be taken easily, yielding a power 
series of £ for the integrand. Therefore, the integration 
over £ is immediately performed. Since all higher order 
in ((cos(<£>))) e f f terms of the expansion vanish at T = T c , 
we get the equation for T c 
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where r\ is the coordination number. The phase correla- 
tor in Eq.(6) is calculated in the clean limit of the ID 
free energy functional, obtained from Eq.(l) by setting 
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0, which returns (see Ref.[22]) 



(cos(p(0)) cos(p(*))) - exp{-|z|/r c }, (7) 
where r c — ft 2 N^iT) / '2toi|£i|T '. Using dimensionless T c - 
shift t = \frjTpE~Z (yf — tTit) anc ^ disorder parameter 
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FIG. 1: The physical solution t(q), giving the T C (W)- 
dependence, within the full range from clean limit (CL: q — 0) 
to the dirty limit (DL: q c = 2/3v3)) is highlighted as the bold 
(blue) curve. Formal solutions of the cubic Eq.(8) are shown 
for completeness. T c (q) vanishes abruptly at q — q c . 
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Eqs.(6,7) yield 



1 = t 2 {l - qt). 
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Expanding the physical solution of this cubic equation in 
the weak disorder regime (small q) the T c -shift obeys 
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(9) 



showing that T c decreases with increasing randomness 
like W 2 . For a pure system Eq.(9) gives the dependence 

1/2 

T c ~ , in agreement with Efetov and Larkin in Ref. 
[18]. This expression shows that even a small interchain- 
coupling sets up an ODLRO in the system, and conse- 
quently, the critical temperature increases with E±. On 
the other hand, disorder reduces T c due to "melting" of 
the order parameter phase coherence between neighbor- 
ing chains. The exact solution of Eq.(8) gives three roots, 
among which the physical one is confined to the finite q- 
range as shown in Fig.l by the bold line. According to 
this solution, the critical temperature decreases mono- 
tonically with increasing q in the interval of < q < 



The SC phase hence suppressed beyond the 

4B i T< 1) 



critical disorder-value W 2 



, being trans- 



formed into a normal metallic phase for W 2 > W 2 . The 
critical temperature drops to zero at W 2 = W 2 with a 



1 
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Thus the 



jump of size AT C = T* - ^ nepE± . ., 

SC-normal metal phase transition appears as a lst-order 
transition. 

In order to find the behavior of t near the disor- 
der limit (DL in Fig.l) {t*,q*} = {V3,-^}, we ex- 
pand St = t* — t in terms of Sq = q* — q, which 
gives St = 3 3 / 4 V&z. In other words, ST C = T c - T* 



behaves as 5T C = 



in the vicinity of the breakdown point {T*,W 2 } — 
l 



4E I T (1) 



A 

3))< i 



Quantum phase fluctuations regime. We shall improve 
the calculation of the phase-correlators by taking into 
account the transverse rigidity of the system, which pro- 
vides a more realistic determination of the transition tem- 
perature in the quantum fluctuation regime. We start 
from the Lagrangian, for simplicity at B = 



dz[hLp s (z)} 2 - F eff {ip} 



(10) 



where ip denotes the time derivative of the phase. 
The dynamical term in C can be interpreted as the 
electrostatic energy of charged wires [19, 21] E e i = 
dz J dz'Cij(z — z')Vi(z)Vj(z') produced accord- 
ing to the first Josephson equation tp = (2e/h)V; 
Cij(z — z') are the specific coefficients of electrostatic 
induction. Rewriting E e i in terms of time-derivative of 
phases, the Fourier transform K(qj_,q z ) of the new co- 
efficients Kij(z — z') = -^C\_^z — z'), has the phys- 
ical meaning of a compressibility. In Eq.(10) we ne- 
glect a spatial dispersion of the compressibility and take 
K(q±,q z ) = K = const. F e ff{ip} is the functional F 
in Eq.(4), written at the saddle point Q yS of the aver- 
aged free energy T . The saddle point of T is found 
to be Cj, g = j2t-Jdz(([l - cos(</?j(z) - ip j+s (z))]}) , 

where T eff = -Tin J Hj, s T> W e ~ F ° ff/T ■ So > Feff has 
to be calculated self-consistently. Note that the model 
would be calculated more rigorously by replacing the 
ID wire with discrete analogue of Josephson-coupled 
cells, and considering a strongly anisotropic 3D Joseph- 
son network. After averaging over disorder, one can 
introduce the 'order parameters' as Ca — {{e llpi )) and 

lab = (( etVi )( e±liPi '))> where ( a and q a ,b are the order 
parameters corresponding to the SC and glassy phases 
with a, b being the replica indices. The model can be 
mapped to the solvable Sherrington-Kirkpatrick model 
for the long-ranged phase-phase correlations, neverthe- 
less a solution of the model for the short-ranged (nearest- 
neighbor) phase-phase correlation case, realized in our 
model, is hard task. 

The Hamiltonian, expressed in terms of the phases 
tfj and canonical conjugate momenta ib as % = 
^2- h J Hjipjdz — C, becomes 
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where IIj = = \hK^(Q)y h and S ct as given by 
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Ej_T eff 



{([1 - cos(.<Pi(z) - <p i+g (z))])) 



(12) 

represents either the dimensionless anisotropy-parameter 
or the transverse rigidity of the random system, while <5o 
in Eq.(12) being the transverse rigidity of the pure system 
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So is a natural small parameter of the quasi-lD supercon- 
ductors, which ensures small interchain-coupling energies 
in comparison with the intrachain Cooper-pair energy. 
Indeed, expressing E±_ through the interchain tunneling 
integral Jj_ [18] as Ej_ ~ J\je F yields S = Jj_/T c (1) . 

The quantum description is realized by expressing y> q 
and Ilq as a linear superposition of Bose operators fe q as 

<p« = (m) 1/2 ^ + b -J and n q = J (S) 1/2 (^- q - 

6 q ). In the framework of the self-consistent harmonic 
approximation (SCHA), we rewrite the expansion of the 
cosine operator in Eq.(ll) in terms of the bosonic particle 
number operator A q = 6 q 6 q [19, 21] as 



7} = ^Mq,T) (6^ + 1/2), 



(14) 



where the eigenfrcquency of oscillation w(q, T) is given 



W (q,T) = oj[q 2 z + S 2 cl e-^'^ T h(2 cosq x cosg^] 1 ^. 

(15) 

We express the amplitude of the frequency to = 
(N^\T)/muKrtyi 2 as Q = 2^aT c (1) r 1 / 2 /? i , where 



(16) 



The parameter a, which is an essential parameter of the 
theory, can assume values between zero and one [18]. 
The factor exp{— S a (g, T)} in Eq.(15) is obtained as 



C /„ T , 2irauj ^ l-cos(q_Lg) / 1 



(17) 



where N is the number of unit cells per volume, and A q = 
{exp (fiw(q, T)/T) — 1} _1 is Planck's distribution func- 
tion for phonons. A physical meaning of exp{— S a (g, T)} 
is an average of cos(<£j(z) — tp$+ g (z)) over all one-phonon 
states, exp{-S Q (g,T)} = ((cos(v?j(z) - <£j+ g (z)))) (T) , 
[19]. Application of the SCHA results in a renormaliza- 
tion of the parameter 5 c i, changing thus the oscillation 
frequency w(q, T) by means of the phase-phase correlator 



5 c 2 ,exp{-^(g,T)}. 



(18) 



Note that the SCHA is valid under the condition 

EJ^IW - ifEq ^g^ W < i, which 

means that few phonons are excited in the system. 

It is easier to sec that F e ff can be calculated as 
T a ff = —ThiTr{e~ n / T } by neglecting the dynamical 
term (K — or a = 0) in Eq.(ll), which gives T e ff = T. 

Let us start with the T = limit: expressing the 
phase-phase correlator e~ Sa ( g, °) in terms of S a (g, 0) = 

te, s ives e " Sa(s,o) = (v(o))° = c> 

which implies that even small interchain-coupling sta- 
bilizes ODLRO, hence also a finite T phase transition 
should exist. In order to get an explicit expression for 
the dependence of S qu on So and on disorder, we have 
to solve the equation S 2 U = 5 2 l e~ Sa ^ s ' ^ together with 
Eq.(12) for S c i, the latter of which also depends on S qu . 
Thus the equation for the reduced transverse rigidity 

3qu — S qu / 5 q u, where Squ = S^ a is the renormalized 
transverse rigidity for the clean system at T = 0, as- 
sumes the form 

{5* qu f- 2a = {S* qu f- a -1 q u, (19) 

where the quantum parameter of randomness q qu reads 

CW 2 ^ 



Qqu 



2E\ 



(20) 



The numerical solution of Eq.(19) is depicted in Fig. 2. 
The reduced T = transverse rigidity S qu {q qu )\T=o is 
shown to decrease with increasing disorder for (fixed) 
a < 1 , and suddenly drops to zero at the critical disorder 
strength q qu = q c qu - Hence the quantum critical behav- 
ior corresponds to a first order phase transition. Fig. 2 
shows how the breakdown point shifts with increasing 
a to larger disorder, and the jump vanishes for a — > 1. 
Eq.(19) becomes linear for a = 1 and gives, by inferring 
the q qu (H / )-relation from Eq.(20), 



Squ(W)\ T =0,a=l = 5 2 



1 - 



CW 2 e F 



2E ± T C 



(1)2 



(21) 



Here, the transverse rigidity 5 qu {W)\T=o,a=i decreases 
linearly with increasing W 2 and vanishes at W ( 



2 _ 



-. The quantum critical behavior in the model 
is however controlled by two parameters, the strength of 
randomness q qu (W) and the parameter of quantum dy- 
namics a. For a < 1, the superconductor- normal metal 
phase transition at T = is always discontinuous, and 
only turns into second-order at a = 1. 

Let us now study the finite T behavior of the trans- 
verse rigidity. The phase transition in a quasi-lD su- 
perconductor occurs at some temperature T = T c when 
the transverse rigidity in the ensemble of phases {<p$(z)} 
vanishes. The energy spectrum co(q±,q z ) of the collec- 
tive excitations is reorganized and the transverse qj_- 
dependent part of w(qj_,q z ) vanishes at T = T c , i.e. 
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FIG. 2: The dependence of the reduced T — transverse 
rigidity S qu (q qu ) on the disorder-strength parameter q qu is 
shown for < a < 1 in steps of Aa = 0.1. At q qu = q qu , 
5 qu (q qu ) drops to zero for a < 1 and vanishes continuously 
only at a = 1. Inserts show the a-variation of the jump 
(upper right corner) and of its position q qu (a) (lower left). 



symmetry breaking occurs in the bosonic excitation at 
T = T c . Inserting the solution of Eq.(17) for T < aT c (1) 
into 5 2 qu {T) = e - s =(s. T ) and using Eq.(18), we obtain 



transition temperature T = T Cl where the phase coher- 
ence between neighboring wires is destroyed. Random- 
ness in the Josephson coupling shifts T c to lower tempera- 
tures and, therefore, the magnetic field parallel to the SC 
wires penetrates easier into the organic superconductor. 
On the other hand, the randomness does not break the 
Cooper pairs, keeping the penetration of a perpendicular 
magnetic field into the SC wires unchanged. 

In this Letter we studied disorder-effects on T c and on 
the diamagnetism of Josephson-coupled quasi-lD super- 
conductors. Interplay of disorder with quantum phase 
fluctuations plays a central role for superconductor- 
normal metal phase transitions in quasi-lD superconduc- 
tors. The quantum criticality is controlled by two quanti- 
ties, namely disorder strength and dynamical parameter 
of phase fluctuations. The present model's quantum crit- 
icality signals the existence of a quantum critical phase 
between SC- and normal phase. Its nature, whether it is 
a "mixing" of a glassy and CDW or SDW phases, needs 
further investigation. 

We thank the DFG for support under grant Op28/7-l. 
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T 6 qu {0) 

T C Squ (T) 



where a new temperature scale is introduced by means 
of T c0 = 5 qu (0)Tc ■ , and C is a constant C ~ 1. In 

exp{— 0/y}, which has a 



terms of, = (s^ 3 
Eq.(22) assumes the form y - 

non-zero solution only for 9 < e _1 . The finite solution of 
this equation vanishes discontinuously at 6 — 9 C = e _1 , 
giving the following value for T c 



T c = T cQ a~^(2/eC) — 



(23) 



The magnitude of the jump in y(0 c ) is e _1 , and hence 
the phase transition is of first-order. 

Meissner effect. The current density is calculated ac- 
cording to -J(2,j) = — T-A-(ln Z(A)). For simplicity we 
present here only the diamagnetic contribution to the i- 
th (i = ||, _L) component of the current 



4ttA; 



where the longitudinal- Ay and the transverse A_l com- 
ponents of the penetration depth are obtained as 



h 2 



((cos(v?j - </>j H 



*)»■ 
(25) 

Although Ay(T) diverges at T = T c (1) due to pair break- 
ing in the SC wires, X±(T) diverges at the global SC 
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